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Abstract 

We consider a non relativistic particle on the surface of a semi-infinite cylinder of circum- 
ference L submitted to a perpendicular magnetic field of strength B and to the potential 
of impurities of maximal amplitude w. This model is of importance in the context of 
the integer quantum Hall effect. In the regime of strong magnetic field or weak disorder 
B » to it is known that there are chiral edge states, which are localised within a few 
magnetic lengths close to, and extended along the boundary of the cylinder, and whose 
energy levels lie in the gaps of the bulk system. These energy levels have a spectral flow, 
uniform in L, as a function of a magnetic flux which threads the cylinder along its axis. 
Through a detailed study of this spectral flow we prove that the spacing between two 
consecutive levels of edge states is bounded below by 2iTaL~ 1 with a > 0, independent 
of L, and of the configuration of impurities. This implies that the level repulsion of the 
chiral edge states is much stronger than that of extended states in the usual Anderson 
model and their statistics cannot obey one of the Gaussian ensembles. Our analysis uses 
the notion of relative index between two projections and indicates that the level repulsion 
is connected to topological aspects of quantum Hall systems. 
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1. INTRODUCTION AND RESULTS 

Recently there has been mathematical progress concerning the spectral properties of dis- 
ordered quantum Hall systems with boundaries. In the theory of the integer quantum Hall 
effect one considers non-interacting electrons confined on the surface of a finite cylinder [1] 
or on a corbino disk [2], submitted to a perpendicular uniform magnetic field of strength 
B and to the potential of impurities of maximal amplitude w. In a classic paper on the 
subject [2] Halperin argued that, at least for strong magnetic field and weak disorder 
(B » w in appropriate units), there exist quantum mechanical states localised near and 
extended along the boundaries of the sample. These states carry a diamagnetic current 
contributing to the total Hall current. Halperin's analysis applies to energies that lie in the 
gaps separating the Landau bands of the bulk disordered hamiltonian, i.e the hamiltonian 
of an infinite two dimensional planar system (with no boundaries). Here we will call this 
part of the spectrum the "pure edge spectrum". Progress towards the characterisation 
of the nature of the pure edge spectrum has been made in recent works for systems with 
one smooth boundary [3], [4], [5]. In the present contribution we obtain new results for 
such systems, which are used in separate work on more realistic geometries involving two 
boundaries [6]. 

We consider the Hamiltonian of a particle on a cylinder of radius ^ thread by a flux 
line with flux $ 

#(<*>) = \vl + \{Py - Bx + |) 2 + W(x) + V(x, y) (1.1) 

where x G R, -| < i; < |, with periodic boundary conditions in the y direction 
^(x,— 77) = \l/(x, I 1 ). The particle is confined to the left half of the cylinder because 
of the external potential W which models the boundary of a " semi-infinite cylinder" . We 
assume that it is continuous, and W(x) = for x < 0, W'(x) > for x > 0, W(x) — > +00, 
x — > +00. For technical reasons we assume a growth of W that is not too fast: we suppose 
that for x > 0, uix 1 < W(x) < u^x 1 , for some < u\ < U2 and 7 > 2. The potential 
of impurities V is piecewise continuous and bounded < w with < w < ^. We 

also suppose that V(x,y) = for x > 0, however our methods can be adapted to a more 
general model where the impurity potential extends inside the region of the boundary. 

We will also use two other Hamiltonians: the "edge hamiltonian" H e (&) obtained from 
(1.1) by removing V and "the bulk hamiltonian" -ff&($) obtained from (1.1) by removing 
W. 

The "semi-infinite planar" case corresponds to L = +00. In this limit the correspond- 
ing Hamiltonians become independent of $ and we denote them H^, H e oo , Hb j00 . It is 
easy to see that H b)00 has gaps G n d] (n + \)B + w, (n + |)S — w[, n E N. A basic fact is 
that for weak enough disorder the "pure edge spectrum" a(H 00 )nG n , n E N is continuous. 
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This result is also proven for W replaced by a Dirichlet boundary condition at x = and 
for smooth curved open boundaries (see [3,4,5]). 

When L is finite G n contains only discrete isolated eigenvalues. We formulate this 
result and all the subsequent ones in the special case n = 0. 

Lemma 1. Let B > 2w. For any < e < ^ - w the set cr(H($)) n G , G =]§ + w + 
e, — w — e[ contains only a finite number of isolated eigenvalues of finite multiplicity. We 
label the eigenvalues of H(0) in Go as -E'i(O) < -E^O) < ... < -EW(O) for some finite N. Any 
Ek(0) G Go can be continued into one or several analytic branches for $ e [0, 

for some small enough <&k > 0. 

The discreteness of the spectrum in the specified interval is non trivial even if the 
circumference of the cylinder is finite because the impurity potential can extend to infinity 
in the direction x — > — oo where there is no confinement. In fact one can see that the 
rest of the spectrum may have dense parts. For example if V is a typical realisation of a 
random potential the Landau bands [(n + |)S — w, (n + \)B + w] have dense spectrum. 
Now let0<5< -f — w — e and A =}B — 5,B + 8[. For L large enough, as long as an 
eigenvalue Ek(<&) € A for some $, then we are assured that it can be continued into an 
analytic branch for the whole interval [0, 27r]. This comes from the fact (see inequality 
(3.15)) that the maximal variation of Ek(&) is 2tt\/3BL~ 1 so that it stays in Go and never 
merges in the Landau bands. 

In the rest of this work we fix e small and < 5 < ^ — w — e, and look only at 
eigenvalues Ek(<&) £ A. Note that as $ varies from to 2n some of the branches may 
move in or out of A. A reformulation of the analysis in [3,4,5] shows that there exists a 
spectral flow which is uniform in L. This is expressed by the following Lemma. 

Lemma 2. Let B > 2w. There exists 5 : wq small enough, L large enough such that for 
w < wo, L > L all eigenvalues E k (&) e A satisfy 



where a is strictly positive independent of L and k, and depends only on W, B, w and S. 

The existence of a spectral flow is equivalent to the presence of a chiral diamagnetic 
current. Indeed by the Feynman-Hellman theorem 




(1.2) 



d 



E k (*)=j k (Q) 



(1.3) 



where 



j fc ($) = i < * fc (^)|(p v - Bx + |)tt fc ($) > 



(1.4) 
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is the diamagnetic current (or edge current) associated to the eigenstate \^ k (&) > corre- 
sponding to the level E k (&). 

The hamiltonians H(<&) and H(<& + 2tt) are unitarily equivalent, the unitary operator 
being multiplication by exp(27ri^-). Thus for each -£&($) which does not merge in the 
Landau bands there must exist some k' such that E k {2%) = E k >(0). From Lemma 2 it is 
clear that k! > k, but this does not characterise completely the spectral flow. Our main 
new result states that k' = k + 1 and characterises the level spacing for the pure edge 
spectrum. 

Theorem 1. Let B > 2w. There exist 5, w$ small enough, Lq large enough such that for 
w < u> , L > L , the branches E k (&) belonging to A for all $ e [0, 2n] satisfy 

E k (2n) = E k+1 (0) (1.5) 

Moreover the level spacing in A satisfies 

^<\E k+1 (0)-E k (0)\< 2 -^- (1.6) 



For the constant a in Lemma 2 and theorem 1 we can take the right hand side of 
(2.29). The important point is that in the lower bound of (1.6) a does not depend on the 
detailled configuration of the impurity potential but only on its maximal amplitude. So 
for a random potential the level spacing is random but our lower bound is non random. 

For the usual Anderson model it is proven that the level spacing of localised states 
satisfies Poisson statistics [7] , [8] and it is numerically established that extended states have 
a level repulsion satisfying the Wigner surmise [9] . Here we have a different situation: the 
states are extended, chiral and have a much stronger level repulsion which makes the level 
spacing very rigid. Let p(E) denote the average density of edge states. We expect from (1.6) 
that, in the limit L — > oo, the rescaled level spacing s = Lp(E k )\E k+ i — E k \ has a histogram 
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p(s) which is a certain broadening of S(s — 1) with a finite support of O(jrj). The level 
statistics cannot follow the Gaussian ensembles and it would be worthwhile to investigate 
this question numericaly for an analogous model on a lattice. It is apparent from the proof 
of theorem 1 that the rigidity of the edge spectrum is related to the topological invariants 
of the quantum Hall effect. Also if the spectral flow would satisfy E k (2n) = E k+n (0) with 
n > 2, it would not be forbidden to have n consecutive levels arbitrarily close. 

We wish to point out that all these features can be checked immediately for a simple 
toy Hamiltonian. Consider a one dimensional chiral particle on a circle of circumference L 
thread by a flux $ 

h(*) = (-id v + *) + v(y) (1.7) 
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The exact spectrum is 



LLP 



L 



2ixm $ 1 f 2 7 
t (d>) = + + dyv(y) (1.8) 



which satisfies (1.2), (1.5), (1.6) and has p(s) = 5(s — l). It is expected that (1.7) is a good 
approximation of (1.1) for distances to the boundary of the order of the magnetic length 

Finaly we recall how it follows from (1.5) that the "edge conductance" of the semi- 
infinite system is quantized (see [1], [2], [4] for similar discussions). Let Pa(&) be the 
projector of -ff($) on an energy range A. The edge conductance may be defined as the 
total edge current per unit energy, 

ae = z5L TaIl 1 ^ " Bx ~ t) Pa ^ ] (L9) 

We assume that for a suitable class of potentials V this limit exists and is independent of $ 
(the flux has no effect for the semi-infinite plane). We expect this assumption to be true for 
typical realisations of random potentials that are ergodic with respect to the translations 
along y. In this case the limit should be equal to ^Av J dx < x,0\(p y — Bx)P< x ,a\x, > 
where Av is the average over the disorder and -Poo, a the projector of onto A. The 
limit of the later quantity when A — > fx has been shown to be an integer if fi is a point in 
the gap Go, by non-commutative geometry techniques applied to the lattice case [10]. In 
the present situation it is easy to see that for A in the first gap of the bulk Hamiltonian 

Hb,oo 

1 \\(Py -Bx- f )Pa(<&)|]i < T^rrWiPy -Bx- ?)Pa(<&)||.||Pa(<&)||i 



\A\L nyry i ~ |A|L L 

< j^|sup| WNl (< V|Pa(*)(# (*) - V)P±(*)\if, >) 1 / 2 TYP A ($) (1.10) 

< j^TrP&W = 0(1) 

Here ||.||i and |.| are the trace and operator norms respectively and we used||AB||i < 
||A||.||5||i for A bounded and B trace class. In the last equality we used that there are 
O(L) states in A because of (1.6) so that the final bound is uniform with respect to L. 
Since we have assumed that a e is independent of by averaging over $ we get 

1 f 2ir d$ dE k (<$>) 1 fc ^ [ 2lT d$ dE k m 

a e = hm — -— / — — > — — = hm — -— > 



l^oo A Jr. 2% ^ d§ l^oo A ^— ' L 2tv d§ 
1 1 JU E k (<i>)eA 1 1 k min JU 

lim — i-j-r V (E k+l (0) - E k (0)) = lim — m(Ek - E k .) = — 
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For the first equality we use (1.3), (1.4) and dominated convergence. To obtain the second 
equality we consider separately the contributions of the eigenvalues with k m i n < k < k max 
such that E k (&) G A for all $ G [0, 2n], and of a finite number of eigenvalues with k < k min 
(resp. k > k max ) which enter (resp. leave) A as $ varies from to 2ty. From (1.6) and 
(3.15) this later contribution is 0(L _1 ). Finaly (1.5) is used in the third equality. Here 
the units are such that e = fi = 1 so tt- = 4-. 

2-k h 

Section 2 contains the proofs of Lemmas 1 and 2 and a third Lemma that is needed for 
the proof of theorem 1 in section 3. The appendices A and B contain technical estimates. 



2. DISCRETENESS OF EDGE SPECTRUM AND SPECTRAL FLOW 
Proof of Lemma 1 

Let D > to be chosen later (large) and Vd(x, y) = V(x, y) for x < —D, Vd(x, y) = 
for x > —D. Then V(x, y) — Vd(x, y) has compact support and a standard argument using 
the resolvent identity implies that the essential spectra of 

H D ($) = H e ($) + V D (x,y) (2.1) 

and 

H(Q) = H D (<S>) + V(x, y) - V D (x, y) (2.2) 

coincide [11]. Therefore if we show that <j{Hd{Q)) nGo contains only isolated eigenvalues 
of finite multiplicity, the same is true for H(&). This will be achieved below using a 
decoupling scheme [12], [13] which proves that <j(Hd(&)) H Go is a small perturbation of 
a(H e (G>)) n G Q . The set a(H e (G>)) consists of non degenerate energy levels e n (^j^ + ^), 
n£N the Landau index and m G Z, where e n (k), fceR the wavenumber conjugate to y, 
are the spectral branches of H e oo . These spectral branches are monotone increasing entire 
functions of k with e n (k) — > +oo for k — > +oo and e n (k) — > (n + \)B for k — > — oo (see for 
example [3]). 

In order to set up the decoupling scheme we introduce the characteristic functions 
Xe(x) of — y < x < +oo and Xb( x ) of — oo < x < — y. Note that Xe( x ) + Xb( ;r ) = 1 f° r 
all x. We also need the monotone and twice differentiable functions J e (x), Jb{x) such that 
J e (x) = for -oo < x < -^P - 1 and J e (x) = 1 for + 1 < x < oo; J b (x) = 1 for 

-oo < x < - 1, J b (x) = for - T + 1 < x < oo. 

We introduce the Green functions G a {z) = (H a ($>) — z)~ l for a = e,b,D and z G C 
in the resolvent set of the corresponding hamiltonian. Since 



H D (<S>)J a = H a ($) J a for a = e, b 
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(2.4) 



following [13] we have 



(H D ($) - z)(J e G e (z) Xe + JbG b (z) Xb ) 

= (H e ($) - z)J e G e (z) Xe + (H b ($) - z)J b G b (z) Xb 

1 1 (2-5) 

= J e X e + J b Xb + '[Pi, Je]G e (z)Xe + ~[pl, J b ]G b (z)Xb 

= 1 + K e (z)+K b (z) 
where K a (z) = \ \p 2 x , J a ]G a (z)x a , a = e, b. Thus 

(H D (Q) - z)- 1 = (J e G e (z) X e + J b G b (z) Xb )(l + K e (z) + K^z))' 1 (2.6) 



In Appendix A we prove the following estimates for the operator norms of K e {z) and K b (z) 
f + w < Rez < 



for ^ + w < Rez < — w (in what follows c is a generic positive numerical constant) 



\\K e (z)\\<f^e-^ D2 (2.7) 

\mz)\\<j^-e-^ D (2.8) 
Oo{Z) — cw 

where 5 e (z) = dist(z, a(H e ($>)) and where 5q{z) = min(|;z — |, \z — ^|). We have to 
take w small enough so that the denominator in (2.8) stays positive. Later on we choose 
z appropriately and D large enough so that both terms become smaller than |. Thus 

(H D ($) - z)- 1 = J e G e (z) Xe + JbG h {z)xb + R(z) (2.9) 

where 



\\R(z)\\< \\G e (z)\\ + \\G b (z) 



(1 - \\K e (z)\\ - \\K b (z 



-l 



(2.10) 



Let m G Z be such that £0(^77^ + f;) is an eigenvalue belonging to a(H e ($)) fl Gq. We can 
choose p > small enough independent of m and L such that the circle C m with center 
e o(^7^ + x) and radius £ encloses only one such eigenvalue. By choosing z in a sufficiently 
thin annulus around C m and D large enough, (2.7) and (2.8) can be made smaller than 
cB?L 2 e -cVBD < 1 At the game time from ^2.10) we have 

\\R(z)\\< C -^e-^ D (2.11) 
P 

so that from (2.9) (H D ($>) — z)~ x is well defined for z in a thin annulus surrounding C m . 
Therefore we can compute the spectral projection Pd(to, <E») of Hd(&) for the interval 
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Im =]eo(^P + f ) - L^oi^T 1 + f) + l[ b y Cauchy's formula. Let P e (m,$) be the 
projector of H e ($) corresponding to the level e (^p + ?)• Thanks to (2.9), (2.11) we 
obtain for D large enough 



HPdK $) -P e (m, $)|| < e-cv^g < ! (212 ) 

P 

This estimate implies that <t(-Hd (<£>)) n/ m contains only one eigenvalue of multiplicity equal 
to one. Note that this conclusion holds for all I m C G . Finally since H e ($) and H b (&) 
have no spectrum in (U m / m ) c n Go we deduce from (2.7), (2.8), (2.9) that Hd(&,L) has 
no spectrum in that same set. Therefore g{Hd{$, L))C\Gq consists of isolated eigenvalues 
of multiplicity one. 

It remains to show that an eigenvalue E k (0) e G can be continued into one or several 
analytic branches E k (&) for $ small enough. In the present case it is sufficient to show 
[11] that (py — Bx) is relatively bounded with respect to H(0). For any ip in the domain 
of H(0) and any complex number z with Imz ^ we have 

\\\{ Py -Bx)^\\ 2 <<^\{H{Q)-V)^> 

=< V>|(#(0) - z)-\H(0) - z)\(H(0) -z + z)ifj> - < iP\ViP > 

< IKi^(o) - ^-^i-n^co)^!! 2 + k|-j|^ll 2 + kj 2 ||(^(0) - ^)- 1 ||.|]^|| 2 + «;1|^|] 2 

-^ l|if(0) ^ l|2+(kl + ^ + w)IH|2 

(2.13) 

This concludes the proof of the Lemma. 

Remark: In (2.13) we can take \lmz\ as large as we wish so the size of the interval of 
analyticity is not limited by the relative bound but rather by the fact that the branch 
E k (&) may merge in the Landau bands (outside of Go) where it may not be isolated 
anymore. Inequality (3.15) shows that for L large enough the maximal variation of E k (&) 
is 2it\ / 3BL~ 1 , so that if E k (§) is contained in A for some $ then it is contained in Go 
and it is analytic for all $ G [0, 2ir]. 

Before presenting the formal proof of Lemma 2 we would like to point out that in fact 
(1.2) is closely related to the ideas in [3] and [4]. Using the unitary translation operator 
x — > x + -gjr and the Feynman-Hellman theorem it is easy to see that 

L-^E k ($) =< * k ($)\(W' + d x V)* k ($) > 

where > is the eigenstate with eigenvalue E k (&). Using the methods of [3] or [4] 

one may show that for E k ($>) e A, \^ k (&) > is mainly concentrated near the region where 
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W'(x) is large so that (1.2) holds provided both V, d x V are small enough. Here we follow 
a different method which is closer to the original argument of Halperin [2] in that it uses 
directly the relation (1.4) instead of (2.13). Only the smallness of V is required. 



Proof of Lemma 2 

The eigenstates \u nm (&) > of H e (&) with eigenvalues e n (^ ri + ^) are of the form 

< xy\u nm ($) >= e i2 ^ y h nm (x) (2.14) 
so that < u nm (&)\(p y — Bx— j;)u n > m '(&) >= for m ^ m' and all n, n' . Therefore writing 

|¥ fc (*)>=|*k(*) >+!**(*)> (2.15) 

where 

+oo 

i*° fc (<*>)>= E c°nuomm> (2.i6) 

m= — oo 
+oo 

l*fc(*)>=E E cri«nm(*)> (2-17) 
n>l m= — oo 

we obtain from (1.3), (1.4) 

, +00 -j- 

L— £7 fc ($)= E |cr| 2 <«o m ($)|fe- J Bx- z )« 0m ($)> 

m= — oo 

+ 2Re < *£(<!>) |(p„ - 5* - > (2 - 18) 

+ <vl>i($)|(^-5x-|)^($)> 

First we show that the last two terms on the right hand side of (2.18) are bounded by the 
norm V3B\ |. The Schwartz inequality implies 

| < *° k m(Py-Bx- > | < \\^l(m-\\(Py-Bx + j)^ k m\ 

<v / 2||vl>i($)||«^($)|^ e ($)^($) >f' 2 



<V2\\^i(m[<nm\H e (m° k m> + 

1/2 



(2.19) 



<*i(*)|ff e (*)*i(*)> 

= V2||*i(*)ll(< ¥ fc (*)|ff e (*)**(*) >) 1/2 

< v^||*i(*)||(^(*) + ^) V2 < V3S||*1(*) 
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For the third matrix element on the right hand side of (2.18) the same method leads to an 
identical estimate. From the Feynman-Hellman formula we have 

< UOm(^)|(p» - Bx - -)«0m($) > = L ^, e ^— + j) 

,,271771 <E> (2.20) 

where eo(fc) is the derivative of the lowest monotone increasing spectral branch correspond- 
ing to the hamiltonian H £j00 . From (2.18), (2.19), (2.20) 



L 1 E ^)^ E l^o(^ + y)-2v^||^) 



m= — oo 

>^ F (M) ^ |cr| 2 -2v / 3S||^($) 



(2.21) 



|m— M|<m 



with the Fermi velocity 



u F (M) =mm| m _ M |<^' (— + -) (2.22) 



The integers M and m will be choosen conveniently below. Writting the Schrodinger 
equation in the form, 



oo + 00 / O <T> \ 

E E 4 m Un(^ + I )-E k ^))\u nm ^)>=V(x 1 y)\^ k ^)> (2.23) 

n=0 m=-oo ^ ' 



and taking the norm on both sides 

oo +oo / _ \ 2 



oo +oo / O rf> \ 

E E \cr\\^(^ + j)-E k ^)) <w* (2.24) 

n=0 m=-oo 



Dropping the term n = 0, using (e n (^f + f ) - E k ($)) 2 > (f - 5) 2 for n > 1 and 
.E fc ($) e A we get 



OO o 



E E icri 2 = n^(*)ii 2 <^^ (2.25) 



n>l m= — oo 



From (2.24) one can also derive a lower bound for ^2\m—M\<m l c fc m | 2- Indeed retaining 
only the term n = and using the monotonicity of eo( 2 ^p 1 + f;) we have 

A(M, m) 2 \4 m \ 2 < w2 ( 2 - 26 ) 

\m— M\ >rh 
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where A(M,m) is the smallest of the two numbers \e Q (^f(M ±rh) + *)| - E k (®)\. Now 
we choose any M such that e (^^- + *) G A and since E k (&) G A we can take m such 
that A(M, fh)> § - 25. Thus 

E ^'"I^T^i (2-27) 

\m—M\>fh ^ 2 ' 

Finaly the normalisation condition for > combined with (2.25) and (2.27) imply 

E (2.28) 

\m-M\<rh y 2 * U > 

From (2.21), (2.25) and (2.28) we have 



L±E k (*)>v F {M) 



3B s w 



2 



v F (M) 7 (f-25)2 



(2.29) 



Clearly vp(M) is a strictly positive number which does not depend on V but only on W 
and B. Therefore (2.29) implies the result of the Lemma for w and 5 small enough. 

It will become clear in the next section that the proof of Theorem 1 requires the 
absence of crossings for the branches E k (&) in A. Since we do not know a priori if this 
is true for H(&), an intermediate step is to construct a suitable perturbation of H(&) for 
which the non-crossing property is satisfied. The perturbation that is added here has the 
effect to lift the degeneracy at each crossing in A in a way that (1.2) still holds for the 
perturbed branches. This is the content of the next Lemma. 

Lemma 3. Fix B, w, 5 and L as in Lemma 2. Assume that V(x, y) is such that the 
eigenvalues Ei(0) are not degenerate. One can construct a finite rank perturbation -R($) 
with ||i?($)|| < L~ 10 such that the spectrum of H(&) = H($) + -R($) in A consists of 
non degenerate eigenvalues forming infinitely differentiable spectral branches which do not 
cross and are labeled as Ei($) with Ei(0) = Ei(0). Moreover the new branches satisfy 

L-^EiW > & (2-30) 
where a is strictly positive and independent of L. 



Proof of Lemma 3. 

Let -Pa($) be the eigenprojector of H(Q) onto A. Then we have 

P A ($)tf($)P A ($) = (2-31) 
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Since the branches Ei($>) are analytic and the eigenvalues are not degenerate for $ = the 
possible crossings are necessarily isolated. Indeed if two branches would coincide on a set 
with accumulation points they would concide over the whole interval [0, 2tt] and therefore 
violate the non degeneracy assumption at $ = 0. Therefore we can assume without loss 
of generality that there is at most a finite number of crossings in A. Let us construct the 
perturbation -R($). First consider the set C of pairs of branches which cross in A (note 
that n branches may cross at the same point and contribute as n ( n ~ 1 ^ pairs). Pick one 
pair of branches in C say (ij) and assume Ei(0) < Ej(0). Suppose they cross at points 
where the label /i takes into account the fact that the branches i and j may cross more 
than once, i.e 

E^) = Ej^) (2.32) 

Let A^ ($) be infinitely differentiable test functions centered at <fr^, with a compact support 
of width Pi and maxo<$<27r |A^($)| < Ai. The real numbers 5\ and Ai will be adjusted 
in a suitable way below. Add to the Hamiltonian H(<&) the perturbation 

= XX(*)0**(*) >< *;(*)| + >< (2.33) 

We take (5\ small enough so that the supports of the test functions do not contain $ = 
and do not overlap. In order to diagonalise the new hamiltonian it is sufficient to work in 
the two dimensional subspace of the branches i and j. The spectral branches of the new 
Hamiltonian do not change for k ^ i,j, whereas for k = i,j they become 

Elm = ±(Em + E j {*) - v /(^($)-E J ($)) 2 + Af J ($)^ (2.34) 

and 

Ejm = I (ex*) + + ^/m^-E^+x^)^ (2.35) 

Since the difference 

Ejm - Eim = ^(Em-E.my + K^Y (2.36) 

is always strictly positive the new pair (ij) is non degenerate for all values of $. Moreover 
by choosing Ai small enough we can make sure that we do not introduce more crossings. 
Therefore the perturbed hamiltonian 

H 1 m = Hm+R 1 m (2.37) 

has a new set C\ of pairs of branches which cross, with one element less than C. One can 
construct in the same way a perturbation i?2($) of (2.37) (with 5 2 , A 2 small enough) so 
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that the new Hamiltonian H 2 (&) = -Hi(<&) + i?2(^) has two less pairs of branches which 
cross than H(&). Since there is at most a finite number of such pairs by iterating this 
construction we end up with the Hamiltonian 

#($) = H(Q) + R p(®) = H (®) + R (®) ( 2 - 38 ) 



of the Lemma, where the sum over p contains a finite number of terms. Note that H(0) = 
H(0) so that the labelling of the Lemma holds. The norm of the total perturbation is 

iijwii < E 11^)11 ^J2 x p ( 2 - 39 ) 

V V 

The condition ||i?($)|| < L~ 10 can always be achieved by choosing at each step 

K < JT§T P ( 2 - 4 0) 

and (3 P < ^. 

It remains to check that (2.30) holds. From the formulas (2.34), (2.35) and Lemma 
2, it is easy to check that at the first step of the construction the new branches have new 
derivatives satisfying 

for all At each step of the construction it is possible to choose test functions such that 

maxo<*< 27r |^Ag($)| < (2.42) 
in a way consistent with (2.40). So at the first step (p = 1) 



d i a 1 

^EU*) >- L - Wl (2.43) 



Of course (2.43) is also valid for the spectral branches of Hi(&) that correspond to k e N. 
Therefore it is valid for all eigenvalues of Hi(&). By iterating the construction we see that 
any branch of (2.38) satisfies 



which implies (2.30). 
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3. RELATIVE INDEX AND LEVEL SPACING 

The main goal of this section is to prove Theorem 1. Let us first outline the strategy 
of the proof. Without loss of generality we can suppose that V is such that E k (0) are non 
degenerate. Indeed if this is not the case one may find a sufficiently small perturbation 
u(x,y), \\u\\oo < L~ 10 such that this hypothesis is satisfied for V + u. If (1.5), (1.6) hold 
for V + u then they hold for V because the perturbation of the discrete levels separated 
by O(L) is at most 0(L~ 10 ). From Lemma 2 we know that for Ek(<&) G A there is a non 
trivial spectral flow: the branches are monotone increasing and since H (0) and H(2tt) are 
unitarily equivalent we must have £7^ (27r) = Ek>(0), k' > k. We want to show that in fact 
k' = k + 1. Let Ep be a single "Fermi energy" lying between two consecutive levels of 
both Hamiltonians Pfo(0) and H(0). Define the integers Qp and Qp to be the number of 
branches of the corresponding Hamiltonians which cross Ep as $ varies from to 2%. We 
will show that Qp = Qf = 1. We know from Lemma 3 that the branches of H(<&) do not 
have crossings, and from the proof of Lemma 1 that the same is true for the branches of 
H D (&). This enables us to relate Qp and Qp to the notion of relative index of a pair of 
projections introduced by Avron, Seiler and Simon [14]. Then by using the fact that the 
Fredholm index of an operator does not change under compact perturbations we deduce 
that Qf = Qp. By explicit computation we can check that Qp = 1 and therefore Qf = 1 
which implies that E k {2ix) = E k+ i{Q). Since the branches of H(&) are a small perturbation 
of those of H(&) we deduce (1.5). Estimate (1.6) is then an immediate consequence. 

In order to make the paper selfcontained we give a short summary of the mathematical 
tools used below, as developed in [14]. Let P and Q be orthogonal projections on a 
separable Hilbert space 7i. The pair (P; Q) is called Fredholm if QP viewed as a map 
from P7i to Q7i is a Fredholm operator. The relative index Ind(P; Q) of the pair is the 
usual Fredholm index of T = QP, that is dimKer(T^T) — dimKer(TT^). One proves 
that (P; Q) is a Fredholm pair if and only if 1 and —1 are isolated finitely degenerate 
eigenvalues of P — Q, when they belong to the spectrum. Moreover one has Ind(P, Q) = 
dimKer(P — Q — 1) — dimKer(P — Q + 1). A useful formula (we use it for m = 0) states 
that if (P — Q) 2m+1 is trace class for some integer m then (P; Q) is a Fredholm pair and 
Ind(P; Q) = Tr(P - Q) 2n+1 , for all n > m. A central result on which we rely is that if 
(P; Q) and {Q; R) are Fredholm pairs and either P — Q or Q — R is compact then (P; R) 
is a Fredholm pair and 



Finaly we note that if (P; Q) is Fredholm then so is {UPU^; UQU^) for any unitary U and 
the relative index remains invariant. Also Ind(P; Q) = — Ind(Q; P). 



Ind(P; R) = Ind(P; Q) + Ind(Q; R) 



(3.1) 



Relation between Qf, Q 



p and the relative index of a pair of projections. 
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We fix Ep G A between two consecutive levels of H(0) and H(2tt) (recall that they 
have the same spectrum). Let Pp,o (resp. Pf,2tt) be the projectors of H(0) (resp. H(2tt)) 
onto the energy range ] — oo, E F ]. We also need the projector on levels £^(0) whose spectral 
branch Ek(&) crosses Ep- Namely 

P c f,o= E P(E k (0)) (3.2) 

E k (0)<E F s.tE k (&)crossesE F 

where P(Ek(0)) is the eigenprojector of H(0) corresponding to the discrete level Ek(0). 
Since Ep G A by taking L large enough we are assured that this sum is finite and that the 
branches crossing Ep remain in A for all $ G [0,27r]. 
Setting PpQ = PFfi — Pp q we have 

Q F = TrP F:0 = Tr(P F , - P%%) = Ind(P F , ; P%%) (3.3) 

We introduce a smooth, monotone increasing function of time (p(t), < t < T, (p(0) = 
and <p(T) = 27r, describing the adiabatic switching of a flux quantum through the axis 
of the cylinder. Let Ut be the unitary time evolution associated to the time dependent 
Hamiltonian H((p(t)). From Lemma 3, as t varies the spectral branches in A do not cross, 
and are monotone increasing. So an application of the adiabatic theorem [15] assures that 
UtPp c U t tends to Pf,2it- Thus there exists some large enough T such that for T > T , 
the pair of projections (Pp c ; U t Pf^Ut) satisfies 

ll^,o- 4^,2^11 <1 (3.4) 

Thus it is Fredholm and Ind(P^ c ; U t Pf,2ttUt) = 0. Since Pf,o — Pp° is finite rank we can 
apply (3.1) to get 

g F = Ind(P F)0 ;P£, C o) 

= Ind(P F , ; U T Pp^U T ) + Ind(U T Pp, 2n U T ; P£ c ) (3.5) 
= Ind(P F , ; UpPp^U T ) 

Finaly let U be the multiplication operator by e l ~ y . Since U does not change the boundary 
conditions and WH(0)U = H(2tt) we obtain the formula 

Q F = Ind(P F , ; U T tfP F , Q UU T ) (3.6) 

The same construction for Hd(0) leads to 

Qf = Ind(P^ ; U^tfp£ Q UU$) (3.7) 
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where Pj? is the projector of Hd(0) onto ] — oo, Ep] and is the time evolution associated 
to the Hamiltonian H D ((p(t)). We remark that the identities of this paragraph can be 
checked by explicit computation for the simple toy Hamiltonian (1.7). 

Remark: In [18] a different relative index for an infinite two dimensional system is studied 
and related to the Hall conductivity viewed as a Chern number. It would be interesting 
to investigate the analogous relationship in the present case with a boundary. 

Equality of Qp and Qp. 

Since V — has a finite support (z — ii^O)) -1 ^ — Vd) is a compact operator 
for z not in a(Hp>(0)). Therefore the resolvent identity and Cauchy's formula imply that 
Pf,o — Ppfi is compact. Thus the pair (Pf,o; Pp,o) ls Fredholm and we can apply (3.1) to 
get 

lnd{P F ^U ] T U ] Pf^UUt) 

= Ind(P Fi o; Pj? ) + Ind(Pj? ; Vpfl P Ffi UU T ) 
= Ind(P F)0 ; P£ ) + Ind(P^ ; tf T tfP£ UU T ) 
+ Tnd(U%,Ul Pg fi UU T ; U ] T UP Ffi U^U T ) 

The first and third terms in the last equality of (3.8) cancel. Thus 

Qf = Ind(Pi? ; V$rfP$ tQ UU T ) 

= Ind(P° UU T P° \P° H -> P^H) 

where in the last line we introduced the Fredholm index of Pp UUTPp viewed as a map 
from Pp H to itself (H the Hilbert space of the cylinder). From Dyson's equation 

r T 

pF,oUUtPF, -Pf,oUU$P% = dsP° Q UU°_ s {V-V D )U 8 P° (3.10) 

Jo 

Therefore the Hilbert-Schmidt norm of the left hand side is smaller than 

rp 

[ ds\\P° Q UU°_ s {V-V D )\\ HS (3.11) 
Jo 

which is shown to be finite in Appendix B. Thus the difference (3.10) is compact and the 
two operators have the same Fredholm index 

Ind(P° UU T P° \P° H -> P° H) 

= 1nd(P% UU$Pg t0 \P% H - P^H) 

which is equivalent to Qf = Qp. 
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End of Proof of (1.5) and (1.6). 

From the analysis of section 2 we know that for D large enough (say D = 0{L)) the 

branches of H e (&) and Hd(&) that belong to A lie close to each other within a distance 
Q( e -cVBLy g ince the 

spacing of the branches of H e ($>) is 0(L 1 ) it follows that Qp = 1 
and therefore Qp = 1. Thus E k {2Ti) = E k+1 (0) and since there exists < <F < 2n such 
that 

E k (27r)-E k (0) = 2n^($) (3.13) 
from (2.30) we get the lower bound 

\E k+1 (0) - E k (0)\ > ^ (3.14) 

Because £7(0) = £7(0), this bound shows that the levels of H(0) (or H(2ir)) are spaced by 
0(L _1 ). Using the spectral flow of H(&), together with the facts that the levels of H(<&) 
and H(<&) are separated by 0(L~ 10 ), and that - is strictly positive, one deduces that 
necessarily E k {2-K) = E k +i (0). Then proceeding as in (3.13) and (3.14) we obtain the 
lower bound (1.6). Finaly the upper bound is a consequence of 

dEu - - 

Ll li m =< * k ^\ p y ~ Bx+ L^ k ^ ] > 

<\\y k ($)\\.\\( Py -Bx+®)* k m (3.15) 

< (< tt fc ($)|2tf($)|* fc ($) > - < * fc ($)|2F|* fc ($) >)* 

< (2£ fc ($) +2w)5 < (3S)3 
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APPENDIX A 

We start with a sketch of preliminary estimates for the Green function of the pure magnetic 
problem on the cylinder of circumference L, 

H m = lpl + l(Py-Bx+j) 2 (A.l) 

Using the spectral decomposition of the Green function Gq(z) = (Hq(3>) — z)~ x on a basis 
of eigenfunctions 

c i3¥tt V»,m(x) (A2) 



and the Poisson summation formula we obtain 

+00 



<x,y\G ($)\x',y'>= £ e^W"^) < x , y - m L\G ,M W, V > (^3) 



m= — 00 



where Go, 00 (z) is the Green function of the pure magnetic problem on the infinite two 
dimensional plane. In the Landau gauge (r = (x, y)) 

<v\G^{z)W >=|r(i-|)tf(±-|,l,||r-^ 

(A4) 

The presence of the Euler V function indicates that the Landau levels remain unchanged 
on the cylinder, and U is the Kummer function [16]. By using some technical estimates 
as in [17] one may show that for y < Rez < ^ the absolute value of (A. 3) is bounded 
above by the simple expression 

^Le-fl^'l 2 V S(x-x',y-y' -mL)e-f(y-y'- mL ^ (A5) 
d °( z > m=-l,0,+l 

where c is a numerical constant independent of B and L. The factor S comes from the 
logarithmic divergence at coincident points 

S(x - x',y - y') = 1 for — |r-r'| 2 >l 

B 2 (A6) 

= In — |r — r'l 2 otherwise 
2 1 1 

3 

A bound similar to (A. 5) holds for \d x < r\Go j00 (z)\r' > |, with cB replaced by cBz and 
S replaced by ^-r'^ 2 wnen -f |r — r'| 2 < 1- The important feature for the subsequent 
estimates is that all the above singularities are integrable. In what follows c denotes a 
generic numerical positive constant. 

Estimate of ||X e ||. 
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From the resolvent identity 

Ke{z) = l[p2, J e ]Go(z) X e + \\vl, Je]G (z)WG e (z) X e (A.7) 

Evaluating the commutator, and using |]G e (<z)|| < S e (z)~ 1 we find 

\\K e (z)\\ < l\\J"Go(z)xe\\ + \\J' e d x G (z)xe\\ + ^{z)' 1 (I \J"G (z) W\ \ + \\J' e d x G (z)W\\) 

(A.8) 

Estimate (2.7) follows from the fact that all norms on the right hand side of (A.8) involve 
matrix elements of Gq{z) and d x Go(z) separated by a distance at least equal to ^. We 
use the estimate (A an operator with kernel A(r, r')) 

||A|| < max^sup J dr\A{v, r')|; sup J dr'|A(r,r')|^ (A. 9) 

For the first norm we have 

3D ii L 3D II 

rfx / dy^(x)| < r|G (z)|r' > \ Xe (x') < — — / dxe~^~ x I Xe (x') 

cy/BL BD 2 
~ So(z) 

(A10) 

In the first inequality we used (A. 5) and in the last one we use the fact that \x — x'\ > ^. 
On the other hand 

J'J(x) [ dx' [ 2 dy'\ < r\G (z)\r' > \ Xe (x') 



EL I — k. 

2 2 



(All) 

o {z) J_d d (z) 

Thus \\J'JGo(z)xe\\ < J^z) e ~ cBL>2 ■ For the term involving d x G (z) the estimates are 
similar. The terms involving W lead to the same estimates provided 

/ dxe-f I^'IV^') and J^'(x) / cfe'e"* I^VCx') (A12) 

J-^-i Jo 

are bounded by (9(exp(— cBD 2 )). This is the case for the class of functions W{x) that 
grow polynomialy as x — > +oo. 

Estimate for II iff, II 
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First we sketch the derivation of an estimate for the kernel of G b (z) and its derivative 
for z in the gap of a(H b (<&)). 



< r\G b (z)\r' >=< r|G (*)|r' > + E / dri " / dr - < r l G o(*)|ri > V(n) 



(A13) 



x < v 1 \G (z)\v 2 > V{r 2 )...V(r m ) < r m |G (*)|r' > 



Here the range of the integrals over xi,...,x m is ] — oo,+oo[, and that of yi,---,y m is 
[— -§]• In order to extract the decay for \x — x'\ large from (A. 13) and (A. 5) we use, 
from B\x - x'\ 2 > 2VB\x - x'\ - 1, 

e -f (|x-x 1 | 2 + |x 1 -x 2 | 2 + ... + |x m -x'| 2 ) < e -^(\x-x 1 \ 2 + \x 1 -x 2 \ 2 + ... + \x m -x'\ 2 ) 
x e -^f-(\x-x 1 \ + \x 1 -x 2 \ + ... + \x m -x'\) e ^ (AAA) 

rn_ _VB_\ I\ _B_i 

< e 16 e 8 1 'e 16 



,n p -^-\x-x'\-^(\x-x 1 \ 2 + \x 1 -x 2 \ 2 + ... + \x rn -x'\ 2 ) 



Thanks to (A.5), (A.13), (A.14) we obtain for § \x - x'\ > 1 

| < r|G 6 (,)|r' > | < ^re-f I— 'I 2 + £ 

< ^ e -^l— 'I 

5o(z) — cty 



(A15) 



This bound is valid as long as w is small enough. Clearly from (A.13), following the same 
steps, we obtain a similar inequality, with cB replaced by CB2, for \d x < r\G b (z)\r' > | if 
§\x - x'\ > 1. 

To estimate ||-K&|| we have to compute the norms on the right hand side of 

\\K b \\ < ^\\JSG b (z)xb\\ + \\J' h d x G b {z) Xb \\ (A16) 

This can be done easily using (A. 9), (A. 16) and the bound (A. 15) together with that on 
the derivative. Then one finds 

\\K h \\< f* L e~^ D (AA7) 
do(z) — cw 
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APPENDIX B 



By Cauchy's formula, and the resolvent identity 

Pj? = [ dz l——=f dz lr7^+ [ dz L—(W+V D ) 

'' u J Tf z-H d (0) J Tf z-H o (0) J Tf z-H D (0) y 'z-H o (0) 

(B.l) 

where the contour IV encloses the part of the spectrum of Hd(0) lying below Ep. Setting 
g = UUj?- s (V — Vd) we have for the Hilbert-Schmidt norm 

\\Pf,o9\\hs < |r F |sup z6rjr || ^_^ ff|| g 5 

+ dist(J, r ^k(o))) ( sup - r - 1 1 W ^kw) gl lHS + wsuPzerp 1 1 ^kw) gllHS ) 

(B.2) 

Here is the length of the contour which is finite because the spectrum is bounded 
below. Since V — Vd has compact support, g is a square integrable function on the 
cylinder. Therefore from the bound (A. 5), (A. 6) on the kernel of (z — Hq) -1 it is easily 
seen that all the Hilbert-Schmidt norms in (B.2) are finite. These norms can be bounded 
above uniformly in < s < T, and the supremum over z stays finite as long as the contour 
does not touch a Landau level. Therefore (3.11) is finite. 
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